U(l) spin liquids and valence bond solids in a large- N three-dimensional Heisenberg model 



O 
O 

> 
O 

in 

(N 



Jean-Sebastien Bernier, 1 Ying-Jer Kao, 1 and Yong Baek Kim 1 ' 2 

1 Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 
2 School of Physics, Korea Institute for Advanced Study, Seoul 130-722, Korea 
(Dated: February 2, 2008) 

We study possible quantum ground states of the Sp(iV) generalized Heisenberg model on a cubic lattice 
with nearest-neighbor and next-nearest-neighbor exchange interactions. The phase diagram is obtained in the 
large- N limit and fluctuation effects are considered via appropriate gauge theories. In particular, we find three 
U(l) spin liquid phases with different short-range magnetic correlations. These phases are characterized by 
deconfined gapped spinons, gapped monopoles, and gapless "photons". As N becomes smaller, a confinement 
transition from these phases to valence bond solids (VBS) may occur. This transition is studied by using duality 
and analyzing the resulting theory of monopoles coupled to a non-compact dual gauge field; the condensation 
of the monopoles leads to VBS phases. We determine the resulting VBS phases emerging from two of the three 
spin liquid states. On the other hand, the spin liquid state near J\ w J2 appears to be more stable against 
monopole condensation and could be a promising candidate for a spin liquid state in real systems. 
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I. INTRODUCTION 

The search for quantum spin liquid states in frustrated mag- 
nets has long been of great interest for condensed matter 
physicists. These exotic states exhibit fractionalized quan- 
tum numbers and may play important roles in understand- 
ing properties of novel complex materials such as high-T c 
cuprates. In the spin liquid states, there exist deconfined 
spinons with spin quantum number S — 1/2, as opposed 
to the more familiar magnon excitations in the Neel state, or 
the spin triplet 5 = 1 excitations in the paramagnetic spin- 
Peierls states. 1 Previous studies of frustrated Heisenberg mod- 
els on two dimensional lattices have revealed possible quan- 
tum spin liquid phases corresponding to the deconfined phases 
of Z2 gauge theoriesi^iiii^. These phases are characterized 
by gapped spinons and gapped ^-vortices that carry the Z2 
fl UXi 6 i 7 J s i 9 i io i ii 4 i2 j-) ue i the gapped spectra of both spinons 
and Z2 -vortices, it is difficult to observe an effect on measur- 
able "local" quantities that distinguish the Z2 spin liquid from 
a regular paramagnetic state. Theoretical proposals" have 
been presented to observe the ^-vortices; however, to this 
date, experimental attempts have shown negative results^ 

U(l) analogues of these phases, the so-called U(l) spin liq- 
uids, are not stable due to instanton effects in two dimensions. 
On the other hand, deconfined phase of U(l) gauge theory 
can be stabilized in three dimensions (3D)iiia and the result- 
ing U(l) spin liquid states support gapped spinons, gapped 
monopoles, and gapless "photons" J&ii Experimentally, these 
phases could be more easily observed since the gapless "pho- 
tons" would show up in the low-energy physical properties 
of the system such as thermodynamics and energy/entropy 
transport. For example, an additional T 3 contribution to the 
specific heat from "photons" could be observed. There have 
been recent theoretical studies of three dimensional models 
in which possible U(l) spin liquid phases are identifiedjA&Ai^ 
These models, however, are constructions in the limits where a 
large degeneracy in the ground state manifold is first achieved 
by hand; for example, in Ref [jjj one needs to first take the 
easy-axis limit of the original Heisenberg model to facilitate 



further analysis. One desires to obtain U(l) spin liquid phases 
without taking these special limits. 

In this paper, we consider the Sp(iV) generalized three- 
dimensional frustrated Heisenberg model on a cubic lat- 
tice with nearest-neighbor(NN) J\ and next-nearest-neighbor 
(NNN) J2 interactions, using the bosonic Sp(A r ) generaliza- 
tion of the physical spin SU(2) = Sp(l) symmetry. Gen- 
eralized spin operators can be expressed in terms of boson 
operators b\ a , bi a at each site i where a = 1, . . . , 2N la- 
bels the Sp(N) index, and the constraint n& = b\ a b la is im- 
posed to fix the number of bosons per site (rif, — 2S for 
N = 1 or the SU(2) limit). This model has an advantage 
over other models as the spin-rotational symmetry is retained, 
and there is a direct connection with the microscopic model. 
We study the large- limit of this model and examine the 
finite- fluctuations via gauge theories. The phases obtained 
in the large- limit could also be relevant to physical sys- 
tems whose microscopic Hamiltonians are "near" the parame- 
ter space of the original Hamiltonian. This idea, for example, 
has been used to explain the observed spin-correlations in the 
2D spin-liquid state in Cs 2 CuCl4 i 21 ' 22 i 23 Moreover, recent ad- 
vances of trapped atoms in optical lattices 24 allow one to con- 
struct atomic systems with enhanced symmetries like Sp(iV) 
or SU(iV); 25 large-TV theories would provide useful informa- 
tion about the ground states of such systems. 

At the mean-field level, we solve the Hamiltonian in the 
N — > 00 limit with k — n^/N fixed. In general, smaller 
(larger) k corresponds to more (less) quantum fluctuations. 
As a result, magnetically long-range-ordered (LRO) phases 
appear in the large k limit and quantum-disordered param- 
agnetic phases arise in the small k limit. Moreover, differ- 
ent LRO and paramagnetic phases show up depending on the 
ratio J2/J1. We find three different U(l) spin liquid states 
obtained by quantum-disordering LRO states with the order- 
ing wave-vectors (tt, tt, tt), (0, 7r, tt) and (0, 0, tt). These U(l) 
spin liquid phases are stable in the large- N limit due to gapped 
monopoles and still have short-range-ordered (SRO) spin cor- 
relations at the corresponding wave-vectors. 

As N becomes smaller, however, monopoles may condense 
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by closing the gap, and confined phases such as valence bond 
solids (VBS) may arise. We examine possible monopole con- 
densation patterns and find resulting VBS states that may 
emerge from the U(l) spin liquid phases. Different VBS 
phases are obtained from the U(l) spin liquid phases with the 
SRO correlations at (n, tt, tt) and (0, 0, tt). On the other hand, 
monopole condensation is not found in the U(l) spin liquid 
with (0, tt, tt) SRO in the simplest consideration of monopole 
condensation. This may imply that one needs to go beyond 
the approximation schemes we used in the analysis, or this 
U(l) spin liquid phase near J x w J 2 is more stable against 
monopole condensation. If the latter is the case, this may be 
a more promising candidate for a U(l) spin liquid state in real 
materials. 

The rest of the paper is organized as follows. In Sec. |HJ 
we introduce the formalism used in this work. In Sec. [TO] the 
large-TV mean-field phase diagram is explained. In Sec. IIVI 
fluctuation effects about the mean-field states are considered 
via U(l) gauge theory and its dual theory of monopoles cou- 
pled to a noncompact U(l) gauge field. Here, three U(l) 
spin liquid states are also identified. In Sec.|V] possible VBS 
phases that may arise at finite-TV are obtained using the dual 
theory of monopoles with frustrated hopping on the dual lat- 
tice. We conclude in Sec. I VII 



the constraint to rib = b\ a bi a — 2NS, where a — 1, ...,27V is 
the Sp(TV) index. For the physical case, TV = 1, S takes half- 
integer values. Following these transformations, the Sp(TV) 
Hamiltonian becomes 



H = £ W*fil>\ a bV)(J-'%b jS ), 



(3) 



where J a ^ — J a p = —Jp a is the generalization of the anti- 
symmetric tensor of SU(2); it is a 2TV x 2TV matrix that con- 
tains TV copies of e along its center block diagonal and van- 
ishes elsewhere^ 

In the TV — > oo limit at a fixed boson density per flavor, 
rib/N = 2S = k, we obtain to leading order in 1/TV a mean 
field theory for S = k/2. The fluctuations about the mean 
field solution give rise to a gauge theory. 2 The mean-field 
phase diagram at TV — * oo is shown in Fig.^as a function of 
J2/(Ji + J2) and 1/k. At large values of S, various magnet- 
ically long-range-ordered (LRO) phases appear and are rep- 
resented by the ordering wave-vector q = Q2, (ft)- The 
short-range-ordered (SRO) phases at small values of S corre- 
spond to quantum-disordered phases with short-range equal- 
time spin correlations enhanced at the corresponding wave- 
vectors. 



II. FORMALISM 



III. MEAN FIELD PHASE DIAGRAM 



The Hamiltonian for our model can be written as 



H = J x Si 



Si 



(1) 



where Si are S — 1/2 operators at site i. Here J\ > is 
the antiferromagnetic exchange coupling on the NN links, and 
J2 > on the NNN links. The proper description of a frus- 
trated antiferromagnet requires that all spins transform under 
the same representation of a group, and two spins can combine 
to form a singlet. This requirement leads one to consider the 
generalization of the physical spin SU(2) = Sp(l) symmetry 
to Sp(TV) and study the quantum ground states of the Sp(TV) 
Hamiltonian in the large-TV limits 

As the first step toward Sp(TV) generalization, one rewrites 
the SU(2) spin operators using their bosonic representation, 
S, = ^b\ a a^bip, where a,J3 =T, I labels two possible spin 
states of each boson, 6j a , and the constraint rib = b\ a bi a — 
2S must be imposed at each site. The Heisenberg Hamilto- 
nian, an additive constant aside, is then given by 



(2) 



where = 



Ji is on the NN links, and = J2 is on the 
NNN links. Here e a p is the antisymmetric tensor of SU(2). 
One can then generalize this expression to Sp(TV) by formally 
introducing TV flavors of bosons on each site, and by changing 



The mean-field theory can be obtained by decoupling the 
quartic boson interactions in S using Hubbard-Stratonovitch 
fields Qij = —Qji directed along the lattice links. The effec- 
tive action then contains the terms 



s = / dr £ir N\Qij\ 2 -QijJ^bfbl 



+ c.c. 



i>j 



(4) 

where r is the imaginary time and the ellipses represent stan- 
dard terms which impose the canonical boson commutation 
relations and the constraint. 2 At the saddle point of the action, 
we get 



(5) 



The one-site unit cell of the cubic lattice has nine of these 
Qij fields. For larger values of S = k/2, the dynamics of S 
leads to the condensation of the bf bosons and one obtains a 
nonzero value of 



(b?) = 



(6) 



This corresponds to a magnetically-ordered phase. 

The large-TV limit of S is taken for a fixed value of n = 
ni, /N = 2S and, depending on the ratio Ji / J2 and the value 
of K, the ground state of S at T = can either break the global 
Sp(TV) symmetry and posses magnetic LRO or be Sp(TV) in- 
variant with only SRO. We optimized the ground state energy 
with respect to variations in (Qij) and xf for different val- 
ues of J2/J1 and k. We also found that each saddle point 
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FIG. 1: Large- N phase diagram of the Sp(iV) cubic lattice model 
with first and second-nearest neighbor interactions as a function of 
Ji /(Ji + J 2) and 1 /k. The LRO phases break spin-rotation symme- 
try. The spin order is collinear and commensurate in the (ir, n, tt), 
(0, tt,tt) and (0, 0, tt) LRO phases while it is helical and incom- 
mensurate in the (q, tt, tt) and (0, q, tt) LRO phases. The SRO 
phases preserve spin-rotation invariance. Notice that the (q,TT,n) 
and (0, q, tt) LRO phases do not have SRO counterparts. [Dashed 
line: first order transition. Solid line: continuous transition.] 

may be described by a purely real (Qij ) . The phase diagram 
is summarized in Fig. [0 Notice the transition between dif- 
ferent phases are in general second-order (solid lines), with 
the exception of two first-order transitions (dashed lines). The 
various magnetically ordered and paramagnetic phases are de- 
scribed in details as follows. 



A. Magnetically ordered phases 

The magnetically ordered phases are characterized by the 
finite condensate xf ^ 0. In this model, we find the following 
LRO states. 



1. (tt, tt, tt) LRO state 

This is a long-range-ordered state in which (S^) is 
collinearly polarized in opposite directions on two inter- 
penetrating cubic sublattices (Fig. [2). A gauge can be chosen 
in which the expectation values of link variables are nonzero 
and equal on horizontal and vertical links, while the values on 
diagonal links are zero. 

2. Helical (q, tt, tt) LRO state 

This helically-ordered phase is characterized by nonzero 
values of (Qij) on all NN links and on four faces of the cube 
(each face has two NNN links (Fig. [3})- In the appropriate 
gauge, two of the three possible NN bond directions present 
the same expectation value. Due to the three rotation axes of 




(7r,7r,7r) 

FIG. 2: (tt, tt, tt) phase. All NN bonds have nonzero values. 




FIG. 3: (q,n,n) phase. All NN bonds have nonzero values and 
NNN bonds are nonzero on four of the six possible faces. 



the cube, there are three possible choices for the nonequiv- 
alent NN link direction: (q, tt, tt), (tt, q, tt), (tt, tt, q). This 
choice dictates on which four faces will lie the NNN bonds 
with nonzero expectation value. It is worth mentioning that 
all nonzero diagonal bonds have the same strength. Also, this 
phase has a long-range incommensurate spin order and the 
spin structure factor peaks at the incommensurate wave vec- 
tor (g,7r,7r). 



3. (0, tt, tt) LRO state 

This magnetically ordered phase has nonzero bond expec- 
tation values for two of the three possible NN bond directions. 
Moreover, NNN links are nonzero and equal on four faces of 
the cube (Fig. [4}- Once again, there are three possible choices 
of direction for the zero NN bond since these three configura- 
tions are interchangeable under rotation around the Cartesian 
axes of the cube. 



4. (0, 0, tt) LRO state 

This phase is characterized by nonzero NN bonds in only 
one of the three possible directions while NNN links are 
nonzero on four of the six faces (Fig.|5Jl. There are three pos- 
sible choices for the nonzero NN bond direction. 



4 




(0,7T,*) 

FIG. 4: (0, tt, tt) phase. NN bonds are nonzero in two directions and 
NNN bonds are nonzero on four out of the six possible faces. [Note: 
dashed lines are only a guide to the eye.] 




(0,0,7T) 

FIG. 5: (0, 0, 7r) phase. Only NN bonds in one direction can be 
nonzero while NNN links are nonzero on four of the six faces. [Note: 
dashed lines are only a guide to the eye.] 

5. Helical (0, q, tt) LRO state 

This helically-ordered phase presents nonzero expectation 
values for all NNN bonds, and two of the NN bond directions 
are also nonzero (Fig. [6}. Diagonal bonds do not exhibit the 
same strength, and the NN bond values continuously decrease 
and reach zero at J\ = 0. Moreover, this phase exhibits a 
long-range incommensurate spin order and the spin structure 
factor peaks at the incommensurate wave vector (0, q, tt). 




FIG. 6: (0, q, tt) phase. NN bonds are nonzero in two directions and 
decrease continuously as Ji decreases. All NNN bonds are nonzero. 
[Note: dashed lines are only a guide to the eye.] 



B. Paramagnetic phases 

The paramagnetic phases are characterized by the vanishing 
condensate xf = and a gapped spinon excitation spectrum. 
We find three distinct SRO states. 



1. (tt, tt, tt) SRO state 

This state is obtained by quantum disordering the (tt, tt, tt) 
LRO state. The expectation values of Qij have similar struc- 
tures as those of its ordered counterpart. In this state, all spin 
excitations are gapped. 

2. (0, tt, tt) SRO state 

Once again, this state is obtained by quantum disordering 
the corresponding LRO state. The expectation values of Qij 
have similar structures as those of its ordered counterpart. Al- 
though all spin excitations are gapped in this state,it is not 
symmetric under tt/2 rotation around two of the Cartesian 
axes (y and z). 

3. (0, 0, n) SRO state 

This state is the quantum disordered counterpart of the 
(0, 0, w) LRO phase. The expectation values of the four NN 
links, not strictly null, are very small and will reach zero as J\ 
decreases. Even though all spin excitations are gapped in this 
state, it is not symmetric under tt/2 rotation around two of the 
Cartesian axes (x and y). 

C. Analysis beyond mean-field for the paramagnetic phases 

It is important to note that these three SRO phases possess 
bipartite lattice structure. To get a clearer picture, we label 
"up-spin" sites, where bosons carry a +1 "charge", with A, 
and "down-spin" sites, where bosons carry a —1 "charge", 
with B. Since all non-zero Q-fields are on bonds between A 
and B sublattices, these fields will form "charge"-zero objects. 
In other words, the mean field parameter (Qij) is invariant un- 
der a global U(l) gauge transformation: 

h -► be**, 
bj - bje-^, 

Qij * Qij: (7) 

where i 6 A and j 6 B. Thus, unlike the Z2 spin liquid cases 
in two-dimensions, the U(l) symmetry of the action is intact. 
These phases would exhibit low-energy gapless U(l) gauge 
excitations ("photons"). 

However, to correctly capture the structure of these phases 
beyond the mean-field level, we also need to consider the con- 
tribution from topologically non-trivial gauge-field configura- 
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tions, namely monopoles. In particular, the interference ef- 
fects between monopole events due to Berry phases lead to 
different quantum ground states. In general, two possible 
classes of paramagnetic ground states exist. In the first class 
of paramagnets, none of the symmetries of the Hamiltonian is 
broken. The spins are paired into valence bond singlets which 
strongly resonate between a large number of possible valence 
bond configurations. Such a state is a spin liquid or a res- 
onating valence bond liquid. In the second class of paramag- 
nets, the valence bond singlets spontaneously crystallize into 
a configuration which necessarily breaks a lattice symmetry. 
This configuration is a valence bond solid. We can distinguish 
these two classes of paramagnets using the expectation value 
of the monopole field as our order parameter. Monopoles are 
gapped in a spin liquid, and condense in a valence bond solid. 

Consequently, in the next sections, we investigate under 
which conditions the three paramagnetic phases are spin liq- 
uids or valence bond solids. First, in section IIVI we write 
down the action for such a system and find in which limit we 
get spin liquids. Then, in section [V] we study the effect of 
monopole condensation and its consequence on the structure 
of the SRO phases. 



However, the presence of the background charge complicates 
the situation. The Berry phase term corresponds to a dual 
magnetic flux emanating from the center of each dual cube. 
This flux alternates in sign following a pattern particular to 
each SRO phase. Thus, we expect to get, after a series of dual 
transformations on S, a theory of monopoles with frustrated 
hopping coupled to a non-compact gauge field. In the remain- 
der of this section, we first explicitly show how to obtain the 
aforementioned monopole action. Then we present the condi- 
tions under which these SRO phases are U(l) spin liquids. 



A. Duality transformations: path towards "frustrated" 
monopole action 

As a first step, we replace the cosine interaction by a Villain 
sum over periodic Gaussians such that 



exp (k cos(A p Ajv - A„A^) 



E ex p ( - f( A ^ - A » A ^ - 27n ^) 2 )' < 9 > 



IV. MONOPOLE ACTION, BERRY PHASES, AND U(l) 
SPIN LIQUIDS 

The SRO phases described above have gapped spinons, 
hence we can safely integrate out spinons and study the result- 
ing compact U(l) gauge theory. 26 This analysis was used suc- 
cessfully in (2+l)D to study antiferromagnets on the square 27 , 
Shastry-Sutherland 3 , and checkerboard 4 lattices. Thus, our 
starting point is the action of a compact ?7(1) gauge theory 
in (3+l)D on a space-time lattice in the presence of a Berry 
phase term. Physically, the Berry phase term can be assimi- 
lated to a static background gauge charge of strength ±1 on 
the sites of the spatial cubic lattice. This action is given by 

/■2jt 



Z A = / V[A m ]e^[K V cos(A^ 3 v - A v A jtl ) 
- i2S^2r}jA jT 



(8) 



where j denotes the sites of the space-time lattice, and the 
cosine term represents the conventional Maxwell action for 
a compact U(l) gauge theory. It is the simplest local term 
consistent with the gauge symmetry A, M — * Aj^ — 4>j — <frj+^ 
and it is periodic under Aj^ — > Aj^ + 2ir. The last term is 
the crucial Berry phase leading to large cancellations between 
different paths. The static background gauge charges, r]j, can 
be ±1 depending on the underlying magnetic structure. 

In (2+l)D, the presence of an alternating static background 
charge lead, in confined paramagnetic phases, to broken trans- 
lational symmetry. Here, our aim is to extend this analy- 
sis to (3+l)D paramagnets by performing duality transforma- 
tions on S. Neglecting the Berry phase term, the remaining 
compact U(l) gauge theory becomes a theory of point-like 
monopoles coupled to the dual non-compact U(l) gauge field. 



where j denotes the sites of the dual lattice. We then rewrite 
this expression using Poisson resummation formula 



exp (k cos(A jU A ;7 v — A„Aj 

e J — OO 



/ K \ 

exp ( ~ y ( A ^ A ^ ~ A » A ^ ~ 2lT< h^f + ^fj^^j > 

(10) 

where fj^ lives on the links of the direct lattice. Following 
the change of variable = A M A jv - A v A m - 2-K<fy-^ v , we 
perform Gaussian integrals over the intermediate fields 
and obtain 

exp \ K cosfA^ A jv - A.vAjfM 
f 2 

-E cx p(-^ + */^( A /^> - . (11) 



2K 



Upon performing these transformation on (|8j, the resulting 
partition function is 



[■2TT 

Z A = V[A m ] Eexp( E 



-f 2 

2K 



+ i 



ifjiivXAftAjv — A^Ajfj,) i2S S rjjAj T 



(12) 



where the integer-valued field fj^ v can be interpreted as an 
analogue of the electromagnetic field tensor. For complete- 
ness, it should be noted that we dropped in Eq. (I12> an in- 
significant overall normalization constant. 
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We can then perform on this more amenable expression the 
integrals over A, M , and we find (for S = 1/2) the condition 

This condition can take the form 

A a &j a — IJj 

e abc A b bj c - A T e ja = (14) 

where e JQ is the "electric field" on the links of the direct lattice 
and bj a is the "magnetic field" on the links of the dual lattice. 
Upon solving these two equations, we obtain 



-■j a 



Cja + e abcA b hj c 



(15) 



Thus, A a e° a = rjj, and we can rewrite the partition function 
in terms of ey a and bj a 

j,j,a 

ex P ( " air E (( e °« + e ^bh lc f 



E 



+ (A T h Ja + A a ^f) 



(16) 



The physical properties of the partition function become 
clearer by parameterizing the "frustration" e" a into a curl in 
terms of a new fixed field Yj , . Thus, 



e ja = e abc A b Y lc . 



(17) 

Inserting Eq. (I17> into Eq. (1161 . we can now write the partition 
function in the following form: 

Z A = ^(-^T,((^oA b (h lc + Y Jc )) 2 



+ {A T h la + A a ^Y)y (18) 

To write Eq. Jl 8i in a simpler form, we define a new dual field 

Lq a = hj a + Yj a , and make the gauge choice ipj = —Lq T . As 
a result, our theory is now given by 



Z A = ^ ex P 

(L-, } 

1 3P 1 



0,P<° 



We then promote the integer-valued field L-j p to real values 



with appropriate conditions. To do so, we first use the Poisson 
resummation formula to soften the integer constraint on Lj p . 
The cost of this procedure is the introduction of a new integer- 
valued field J- m \ the monopole four-current. After shifting 



op 



the real field by L 

Z A = 



j ii 



Ljp + Yj p , we obtain 



/CO i 

j,p<" 

^4™\Llo + Y lo )). (20) 



2ni 



~jp ^ op 



Then, to control the fluctuations of ji m \ we add the mass 

OP 

term 



fugacity 



ln(A m ) 



(21) 



3-l> 



Moreover, to explicitly make the the dual theory gauge invari- 
ant and to enforce monopole current conservation, we add a 
second term given by 

V[8t m) } exp (2z^i m) A p ji m) ) , (22) 

JO 1 \ OP J 

where 6j is a U(l) field existing on each site of the dual 

lattice. Note that exp(z6>i"^) corresponds to the monopole 
creation operator and A m is the monopole fugacity. 

Following these two transformations, the partition function 
is now given by 



/OO />27T 
-oo JO 



Ml -S 



(23) 



i4:'} 



5 = ^E ^ l ip A ^>) 2 



0,P<" 



E (^(jj p m) ) 2 



OP 



2i*j£\*df*-Lj p -Yj p )). (24) 



Finally, upon summing over J- m \ and rescaling all fields by 



2-7T, the final action reads 



OP 



Z A = j V[L- Jp ] j V{6^} e- s 



(25) 



o,p<& 



A m ^cos(A P 4 m) -L 7p ~r^ 



(26) 



OP 



It is important to note that, since the static gauge charge at 
the center of each dual cube is ±1, the quantity e abc A b Yj c , 
interpreted as fluxes through the faces of the dual cubes, is 
defined modulo 2ir. 

To summarize our result, the duality transformations on S 
yield a theory of monopoles with frustrated hopping on the 
dual lattice, coupled to a dual non-compact gauge field Lj p . 

The "frustration" is encoded in the field Yj and arises from 
the Berry phase term. 
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B. U(l) spin liquids 

From the final expression of the dual action, we see that, for 
A m close to zero, the monopole field is gapped, and the non- 
compact dual gauge field is free. Thus, the resulting phase has 
gapped spinons, has gapped monopoles, has gapless photons, 
and respects all the lattice symmetries. This phase is a U(l) 
spin liquid. Thus, for small fugacity, all three SRO phases can 
be spin liquids. It is interesting to note that the sole inclu- 
sion of the NNN interaction in the cubic lattice has allowed us 
to find three neighboring spin liquid phases in the parameter 
space. 



V. VALENCE BOND SOLIDS 

When A m becomes of order one, monopole condensation 
gives rise to phases with broken symmetries. In this section, 
we study different valence bond configurations that may arise 
from the three SRO phases. In order to achieve this, we first 
ignore the fluctuations of the dual gauge field Lj p . Once 
the structure of the resulting slow fluctuation description is 
known, we may restore the field Lq . Hence, we begin our 



study with the following action 



S — — A m 



cos 



(m) 



Y- 



3PJ 



(27) 



3P 



We then rewrite the cosine term using exponentials, and pro- 
mote t to a continuous variable. The action now takes the 
form: 



y « + c.c.) 



{RR') 



(28) 



where R labels the sites on the cubic lattice dual to the orig- 
inal lattice of spins. The sites of the dual lattice are at the 
centers of the cubes defined by the planes of the original lat- 
tice, and R = (x, y, z) with x, y and z integers in units of 

the lattice constant (a — 1). Then, assuming that B^/ is a 
slowly-varying function of time, we find 



where $ 



+ c.c. 



\d T 



(m) 
R 



(29) 



is the monopole creation operator. Con- 



sequently, by dropping the constant term from Eq. d29i and 
adding the potential ^ ( | $ fl J 2 ) = r 1 $_r | 2 + /i 1 | 4 + • • • , we 
obtain the continuous-time soft-spin action for the frustrated 
XY model 

S =J dr^l^Rl 2 - ]T (t RR ,$* R $ R ,+c.c.) 



(RR!) 



(30) 




FIG. 7: Flux of ±7r through two of the six plaquettes of the dual 
lattice. There is no flux along the directions for which spins are not 
staggered. 




FIG. 8: A gauge choice for Y^ p that realizes the fluxes shown in 
Fig.0 



where the frustration is encoded in the monopole hopping am- 
plitudes t RR / = te lYltR ' corresponding to the fluxes through 
the dual plaquettes. These flux values depend on the fixed 
field rjj which distinguishes the different SRO phases, and en- 
capsulates the spin staggering in the local collinear order. It 
has the values 




for(0,0,7r) SRO 
for(0,7r,7r) SRO 
for (7r,7r,7r) SRO. 



(31) 



In the remainder of this section, we will find the possible 
monopole condensation patterns associated with the three 
SRO phases. 



A. (0,0,7r)SRO 

To analyze this SRO phase, we first need to find the 
monopole hopping amplitudes. To do so we choose an appro- 
priate gauge (Fig. [8) that realizes the fluxes on the plaquettes, 
as shown in Fig.0 The corresponding monopole hopping am- 
plitude is given by 

te™ v , 



t R ,R+x 
tR,R+y = t, 
tR,R+z = I; 



(32) 



where t and I are real- valued hopping coefficients. These co- 
efficients emphasize the anisotropic character of the monopole 
hopping due to the presence of diagonal bonds on four faces 
of the cube. After diagonalizing the kinetic part of the soft- 
spin action Eq. ( I30> . we find two low-energy modes. The nor- 
malized real-space wave functions associated with these two 



monopole excitations carrying different lattice momenta are 



* 2 (i?) 



(1 + V2) - e^ y 

^4 + 2^2 
((1 + yg) + e^e* 
V4 + 2V2 



(33) 



To characterize the behavior of this system near monopole 
condensation transition, we consider fields which are linear 
combinations of ^1 and ^ 2 . Any such linear combination is 
at the bottom of the monopole band where there is a contin- 
uum of states for the monopoles to condense: 



$(i?) = ai*i(i?) + a 2 ^ 2 (R) 



(34) 



The phase transition in this system can be explored within 
a Ginzburg-Landau theory by treating a\ and a 2 as slowly- 
varying fields. To simplify the next calculations, we define 
two new fields in terms of ot\ and a 2 : 



h = ai 



ia 2 , 



a.\ — ia 2 . 



(35) 



By studying the action of the lattice symmetries, we find that 
the resulting Ginzburg-Landau functional is required to be in- 
variant under the following transformations: 



T, 



R 



I? 



7172,-Rx 



92 
01 

«02 
-i<t>i 

01 

02 

e i7r / 4 0i 
e-™/ A (j) 2 
\ (01 + 
I (01 + 



2 )(-l 
2 )(-l 



) x+z ), 



(36) 

where the tt/2 rotations are about the lattice points on which 
the monopoles reside. These transformations suggest that the 
simplest invariant takes the form 



with 



J(0 1) 2 )=7VX 2 ! 



(37) 



(38) 



where a a are Pauli matrices. We then write down the contin- 
uum action for the two-component complex field 0(i?, r) that 
respects the above symmetries. In doing so, we have restored 
the dual gauge field L p and included some generic kinetic en- 
ergy S L for L p . 

Ssi ow = J drd?R (|(V P - iL p )cb\ 2 + C/(|0| 2 ) + wJ(0)) 

+ (39) 

where t/(|0| 2 ) = r|0| 2 + u\<j)\ 4 + When r < 0, the 
monopoles would condense. From now on, we study confin- 
ing paramagnetic phases originating from the Coulomb phase 




FIG. 9: Picture of a "sheet-like" columnar VBS state obtained 
when w > 0. The monopole density is increased on the x-odd 
planes which is interpreted as having antiferromagnetic bonds be- 
tween NNN spins preferentially crossing these planes. [Note: dashed 
lines are only a guide for the eye.] 



(U(l) spin liquid) by condensing single monopoles. Ground 
states are selected by minimizing wJ(<fi), and the sign of w 
determines the structure of the resulting phase. 

The expectation values N x and N y are sufficient to charac- 
terize each state since the spatial monopole density is given 
by 

\HR)\ 2 - C Q|0| 2 - (1 + V2)((-l)^ + (-lfNy) 

(40) 

where C is a normalization constant. Depending on the sign 
of w, the following phases are possible. 



1. w > 



In this case, there are four ground states 

(JV S) JV„) = (±1,0),(0,±1). 



(41) 



For each state, the monopole density is the same in every other 
plane perpendicular to the x or y lattice axis. For example, the 
state (N x , N y ) — (0, 1) has an increased density on the x-odd 
planes and an decreased density on the x-even planes. In the 
original spin model, these dual planes are crossed by antiferro- 
magnetic bonds between NNN spins, and there is an increased 
bond energy crossing the x-odd planes. Hence, this state cor- 
responds to a "sheet-like" columnar valence bond solid with 
columns of zigzagging dimer chains running along the z di- 
rection. These columns are called "sheet-like" because they 
lie entirely on the yz plane (as illustrated in Fig. [9}. 



2. w < 



In this case, there are four ground states 

(N X ,Ny) = (±1,±1). 



(42) 



For these states, the monopole density oscillates in the x and 
y directions. For example, the state (N x ,N y ) — (1, 1) has 
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FIG. 10: Picture of "cylinder-like" columnar VBS state obtained 
when w < 0. The monopole density is increased on the x and y odd 
planes which is interpreted as having antiferromagnetic bonds be- 
tween NNN spins preferentially crossing these planes. [Note: dashed 
lines are only a guide for the eye.] 



increased density on odd x and y planes and decreased density 
on even x and y planes. In the original spin model, these dual 
planes are crossed by antiferromagnetic bonds between NNN 
spins, and there is an increased bond energy crossing the x 
and y odd planes. Consequently, this state corresponds to a 
"cylinder-like" columnar valence bond solid. As illustrated 
in Fig. El| f° ur connecting sheets of zigzagging dimer chains 
form a cylinder oriented along the z direction. 




FIG. 11: Flux of ±^ through four of the six plaquettes of the dual 
lattice. There is no flux along the direction for which spins are not 
staggered. 



odd ^ 



z even 




FIG. 12: A gauge choice for Y^ p 

Fig.nu 



that leads to the fluxes shown in 



3. Connection to the VBS order parameter 

Modulation of the monopole density in the dual lattice gives 
rise to the modulation of the dual gauge flux in the dual planes. 
The modulation of the dual gauge flux, in turn, corresponds to 
the modulation of static electric field on the links of the direct 
lattice. Following Ref. we can identify the valence bond 
order parameter with these static electric fields. The valence 
bond order parameter in this phase, corresponding to diagonal 
bonds on the xz and yz faces of the cubic lattice, and links 
along the z direction, can be written as 

(-l)V(S r ■ S r + y + i + S r + i ■ S r + y) . (43) 

(-\ys r -s r+ z ) 

Using this observation, we can make the following identifi- 
cation, similar to that in Ref. l20t 

^vbs ~ N = <j^a<j). (44) 

Since the diagonal bonds occur in pairs, the perpendicular 
components of the spin correlations cancel, and we are left 
with the spin correlation along the link directions. Therefore 
the identification of the diagonal bond order parameter with 
the static electric fields in the link may be justified. 

In the above analyses, there is no modulation of the 
monopole density in the z-direction; therefore, there is no va- 
lence bond solid order formation in the z-direction. 



B. (0,7T,7T)SRO 

In this phase, upon choosing an appropriate gauge fFig. 1121 . 
the monopole hopping amplitudes (corresponding to the 
fluxes on the plaquettes as shown in Fig. II 1> are given by 

tn,R+x = -/= (e + le y ) , 

tR,R+y — I, 

t R ,R+z = l\ (45) 

where t and I are real- valued hopping coefficients. After di- 
agonalizing the kinetic part of the soft-spin action Eq. j30l >. 
we find that the lowest energy band is dispersionless in the 
(k x , 0, 0) direction, as shown in Fig.^] This indicates that 
there is a large degeneracy for the lowest lying mode at this 
level of approximation, and no monopole condensation can 
occur. By going beyond quadratic approximation, this mode 
may become dispersive again and we can then determine the 
correct monopole condensation pattern. Similar situations 
have been encountered in the studies of the quantum Ising 
model on a Kagome Iatticep2^2i2i and numerical studies show 
that a valence bond solid with a larger unit cell is possibly 
the ground state. A similar structure might also be the ground 
state, at finite N, for our (0, tt, tt) SRO phase. This situation, 
however, could also be an indication that, at finite JV", the spin 
liquid is more robust in the case of (0, tt, tt) SRO phase than 
it is for the other two SRO phases found in our study. Further 
studies are needed to properly determine the true ground state 
for this phase. 
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E 




kx 

FIG. 13: Band structure given by the kinetic energy of the frustrated 
XY model in the (0, n, n) phase. Notice that the lowest energy band 
is flat in the (k x , 0, 0) direction. The second lowest band occurs for 
the k- vector (k x , n, 0). 




FIG. 14: Flux of ±? through the six plaquettes of the dual lattice. 



Similar monopole hopping problem was studied in the context 
of an interacting boson model. 20 As such, our analysis which 
follows shares some similarities with the discussions in the 
bosonic model, 20 and for completeness, we present relevant 
details of the analysis. 

We diagonalize the kinetic energy and find two low-energy 
modes 



2(3 - a/6) 



1 - (V3- \/2)e l7rz e i7rx -ie l 



2(3 - \/6) 



V2 



-■ (47) 



Any linear combination of these two monopole excitations, as 
Eq. ( I48> . is at the bottom of the monopole band. Hence, there 
is a continuum of states for the monopoles to condense. 



=<PxV 1 (R)+<f> 2 y 2 (R), 



(48) 



Near monopole condensation, we can analyze the transition 
using the Ginzburg-Landau formalism. By studying the ac- 
tion of the lattice symmetries, we notice that the resulting 
Ginzburg-Landau functional is required to be invariant under 
the following transformations: 



T, 



tt /2,Rxy 
ir/2,Rxz 



Y2 

e~ i7r / 4 <; 

<t>*+4>2 

V2 



Y2 
Pi 

01-02 

V2 ■ 



(49) 



Repeating the same procedure, we first find the monopole 
hopping amplitudes (corresponding to the fluxes on the pla- 
quettes as shown in Fig. 1141. Upon choosing the appropriate 
gauge (Fig.EK these amplitudes are given by 



tR,R-i 



tR,R+y = \l-[l — tt 
tR,R+z = 1 



iir(x+y) 



(46) 




zoddff 



z even 
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FIG. 15: A gauge choice for Yj p that realizes the fluxes shown in 
Fig.El 



It appears from these transformations that the simplest invari- 
ant has the form 



y z 



NlNl 



with 



N a {4> u 4>2) = ^a a ^. 



(50) 



(51) 



We then write down the continuum action for the two- 
component complex field (j>(R, t) that respects the above 
symmetries. While doing so, we restore the dual gauge field 
hp and include some generic kinetic energy for this field. The 
resulting action has the same form as the continuum action 
presented for the (0, 0, tt) SRO phase (Eq. (I39», and, as be- 
fore, the ground states are selected by minimizing wK{<f)). 

The expectation values N x , N y and N z are sufficient to 
characterize each state since the spatial monopole density is 
given by 

MR)\ 2 = H 2 + -L[(-iy Nx + (—l) y N y + (~1) Z N Z 

(52) 

Then, according to the sign of w, we can find which valence 
bond solid structures are allowed. 
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FIG. 16: Picture of a columnar VBS state obtained when w > 0. 
The monopole density is increased on the x-even planes which is in- 
terpreted as having antiferromagnetic bonds between NN spins pref- 
erentially crossing these planes. [Note: dashed lines are only a guide 
for the eye.] 




FIG. 17: Picture of a 3D box VBS state obtained when w < 0. The 
monopole density is increased on all x, y and z even planes, and the 
dimers resonate around cubes center where these three planes meet. 
[Note: dashed lines are parts of the real bonds.] 

1. w>0 

In this case, there are six ground states 

(N x ,N y ,N z ) = (±1,0,0), (0,±1,0), (0,0, ±1). (53) 

In a given state, the monopole density is the same in every 
other plane perpendicular to a fixed lattice axis. For exam- 
ple the state (N x ,N y , N z ) — (1, 0, 0) has an increased den- 
sity on the x-even planes and decreased density on the x-odd 
planes, as shown in Fig. In the original spin model, these 
dual planes are crossed by antiferromagnetic bonds between 
NN spins, and there is an increased bond energy crossing the 
x-even planes. Hence, this state corresponds to a columnar 
valence bond solid with dimers oriented in the x direction. 



2. w < 

For this second case, there are eight ground states 

(N x ,N y ,N z ) = -^=(±l,±l,±l) (54) 

where each of the three signs can be chosen independently. 
The monopole density for these states is periodic in all three 



directions and correspond to three-dimensional "box" va- 
lence bond solids. For example, the state (N x , N y , N z ) — 
-^=(1,1,1) has maximal monopole density for even x, y and 
z, as shown in Fig.[l7| In the original spin model, this state 
has dimers resonating around direct lattice cubes surrounding 
the dual lattice points. 



VI. CONCLUSION 

In this paper, we studied the quantum ground states of the 
Sp(7V) antiferromagnetic Heisenberg model on the cubic lat- 
tice with NN (Ji), and NNN ( 3i) exchange interactions. This 
is the simplest 3D model with frustration. We have shown that 
the sole addition of the NNN interaction on the cubic lattice 
gave rise to three neighboring U(l) spin liquids with mag- 
netic short-range correlations at (tt, tt, 7t), (0, tt, tt), (0, 0, tt). 
Moreover, there are indications that (0, tt, 7t) U(l) spin liq- 
uid is more likely to "survive" at finite N values. This model 
has an advantage over other models because of its direct con- 
nection to realistic microscopic models. Thus, verifying the 
results of our model in real systems would be more amenable. 

We first obtained the mean-field phase diagram of this 
model in the N — > oo limit as a function of J2KJ1 + 
J2), which controls the frustration, and 1/k, which con- 
trols quantum fluctuations. For large values of k, we find 
various magnetically ordered phases; (tt, tt, tt), (0, tt, 7t), 
(0,0, tt) collinear-ordered and (q, tt, tt), (0, q, tt) helically- 
ordered states. Upon decreasing k, we find three quantum- 
disordered paramagnetic phases (SRO) with enhanced spin 
correlations at (tt, tt, tt), (0, tt, tt), (0, 0,7r). Since these are 
quantum-disordered states of collinear-ordered magnets, we 
see already at the mean-field level that these phases are possi- 
ble U(l) spin liquids. 

Methods have been developed by Haldane, 32 Read, and 
Sachdev 28 to study (2+l)D quantum antiferromagnets. In this 
work, we have extended these methods to (3+l)D and ana- 
lyzed the effects of singular fluctuations about the mean-field 
states. These monopole events, or the topological defects in 
the U(l) gauge theory, represent tunneling between different 
topologically distinct sectors. We find that, when N — > 00 
(small monopole fugacity), the three paramagnetic phases are 
indeed U(l) spin liquids. At finite N values (large monopole 
fugacity ), the system may become confined due to monopole 
condensation, and we discover various possible VBS ordered 
phases for SRO (tt, tt, tt) and (0, 0, tt) states. However, it ap- 
pears that, when J\ w J2, the corresponding U(l) spin liquid 
with (0,7r, tt) SRO seems to be more stable in the simplest 
consideration of monopole condensation than the two other 
spin liquid states. We would need a more elaborate analysis 
to determine the ultimate fate of this state. One possibility 
would be that it becomes a VBS state with a very large unit 
ce jj i 29 j 30 Another possibility is that it remains a U(l) spin liq- 
uid near J\ w J2 even at finite N. If the latter is the case, this 
U(l) spin liquid may be a better candidate to be observed in 
real systems. One may be able to gain more insights on this 
issue using numerical methods. 

It is worthwhile to mention that the U(l) spin liquids 
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studied previously, in the context of other three-dimensional 
models, correspond to our U(l) spin liquid with SRO 
(tt, 7t, 7t). 18 - 19 - 20 - 33 On the other hand, we find two other possi- 
ble U(l) spin liquids with (0, 7r, tt) and (0, 0, tt) correlations 
within our model. It is also interesting to note that our anal- 
ysis for the quantum-disordered paramagnetic phases at finite 
N shares some similarities with the studies of an interacting 
boson system at half-filling. 20 The nature of the phase transi- 
tions between different U(l) spin liquid states and VBS phases 
is still not clear and awaits further investigation. 
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